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Abstract 

Anyons - particles carrying fractional statistics that in- 
terpolate between bosons and fermions - have been con- 
jectured to exist in low dimensional systems. In the con- 
text of the fractional quantum Hall effect (FQHE), quasi- 
particles made of electrons take the role of anyons whose 
statistical exchange phase is fixed by the filling factor. 
Here we propose an experimental setup to create anyons 
in one-dimensional lattices with fully tuneable exchange 
statistics. In our setup, anyons are created by bosons with 
occupation-dependent hopping amplitudes, which can be 
realized by assisted Raman tunneling. The statistical an- 
gle can thus be controlled in situ by modifying the relative 
phase of external driving fields. This opens the fascinat- 
ing possibility of smoothly transmuting bosons via anyons 
into fermions and of inducing a phase transition by the 
mere control of the particle statistics as a free parameter. 
In particular, we demonstrate how to induce a quantum 
phase transition from a superfluid into an exotic Mott-like 
state where the particle distribution exhibits plateaus at 
fractional densities. 



Introduction 

Usually, every particle in quantum theory is classified as ei- 
ther a boson - a particle joining any number of identical parti- 
cles in a single quantum state - or a fermion, characterized by 
the sole occupancy of its state. The exchange of two fermions 
leads due to the Pauli principle to a phase factor -1 in the total 
wavefunction, while the wavefunction of two bosons remains 
invariant under particle exchange. More than 30 years ago, re- 
searchers proposed a third fundamental category of particles 
living in two-dimensional (2D) systems, "anyons" 1-5 . For two 
anyons, the wavefunction acquires a fractional phase e l ° un- 
der particle exchange, giving rise to fractional statistics, with 
< < 71. For a few years the physics of anyons remained re- 
stricted to the 2D world 6 , until Haldane presented the concept 
of fractional statistics in arbitrary dimensions 7 . 

Anyons in one dimension (ID) are still unexplored to a 
wide extent. Recently, it has been put forward to create frac- 
tional statistics in a ID Hubbard model of fermions with cor- 
related hopping processes 8 . Anyons are realized in this case 
as low-energy elementary excitations. 

In this Article, we introduce an exact mapping between 
anyons and bosons in ID. We show that a Hubbard model of 
anyons is equivalent to a variant of the Bose-Hubbard model 9 
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Figure 1 : Anyon-Boson Mapping and schematic of the proposed 
experiment, (a) Anyons in ID can be mapped onto bosons featuring 
occupation-dependent hopping amplitudes, (b) Assisted Raman tun- 
neling can selectively address hopping processes connecting differ- 
ent occupational states and induce a relative phase, realizing a fully 
tuneable particle exchange statistics angle 6. Energies are not in 
scale. 



where the bosonic hopping amplitudes are state-dependent. 
This conditional-hopping phase factor breaks reflection par- 
ity in the system, which is an important ingredient to real- 
ize fractional statistics 10 . We propose to realize bosons with 
conditional-hopping amplitudes using assisted Raman tunnel- 
ing in an optical lattice 11 ' 12 (OL). We discuss how the di- 
rect control of the statistical phase can induce a quantum 
phase transition from a bosonic superfluid into a Mott-like 
state, exhibiting exotic Mott shells at fractional densities. The 
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"statistical ramp" transmutes bosons smoothly into "pseud- 
ofermions", with anyons as intermediate steps. 

Anyons in ID are defined 13 ' 14 by the generalized commu- 
tation relations 

ajal - e- iea ^-%laj = S jk: aj a k = j e *&<j~ k ) ] a k a h (1) 

where the operators a^aj create or annihilate an any on on 
site y, respectively. The sign function is such that sgn(j — 
k) = for j = k, hence two particles on the same site behave 
as ordinary bosons. In consequence, anyons with statistics 
= 71 are pseudofermions: while being bosons on-site, they 
are fermions off- site. 



Results 

Mapping between Anyons and Bosons 

We introduce an exact mapping between anyons and bosons 
in ID. Let us define the fractional version of a Jordan- Wigner 
transformation, 

( j ~ l \ 
aj = bj exp iQ^tii (2) 

with ni = a\ai = b\b{ the number operator for both parti- 
cle types. Provided that the particles of type b are bosons, 
[bj,bj] = Sjt and [bj,bi\ = 0, we can show that the mapped 
operators a indeed obey the anyonic commutation relations as 
introduced in equation (1). For a proof see Methods. This 
mapping elucidates that anyons in ID are indeed non-local 
quasi-particles, made of bosons with an attached string oper- 
ator. 

Our ultimate goal is to propose a realistic setup for demon- 
strating an interacting gas of anyons in ID OLs. We therefore 
introduce the Anyon-Hubbard model 

L tj L 

H a = - Jj>tfl 7 . +1 +h.c.) + - - 1), (3) 

j Z j 

where J is the tunneling amplitude connection two neighbour- 
ing sites and U is the on-site interaction energy. By inserting 
the Anyon-Boson mapping, equation (2), the Hamiltonian H a 
can be rewritten in terms of bosonic operators: 

H b = -jt(P]b j+ iJ e »J +h.c.) + V - - 1). (4) 

j z j 

The mapped, bosonic Hamiltonian thus describes bosons 
with a occupation-dependent amplitude Je l6n j for hopping 
processes from right to left (j + 1 — >> j). If the target site j 
is unoccupied, the hopping amplitude is simply /. If it is oc- 
cupied by one boson, the amplitude reads Je l ° , and so on. 
The conditional hopping scheme is depicted in Fig. la. We 
emphasize that the non-local mapping between anyons and 
bosons, equation (2), leads luckily to a purely local, and thus 
viable Hamiltonian. 



As expected from anyons, the reflection parity symmetry is 
broken 10 at the level of the commutation relations, equation 

(1) . The fractional Jordan- Wigner transformation, equation 

(2) , transfers this asymmetry also to the bosonic case: The 
resulting Hamiltonian, equation (4), features a phase factor 
acting only on the left site j and thus violates parity. Even in 
the absence of the on-site interaction, U = 0, the exponential 
operator in equation (4) gives rise to many -body interactions, 
as expected for anyons 10 . 

In the limit U/J —> «>, bosons are impenetrable and each 
site contains at most one particle. In this case, the phase fac- 
tor disappears and the bosonic Hamiltonian H b reduces to an 
ordinary Tonks-Girardeau gas 15-18 . However, we consider lo- 
cal occupation numbers beyond the hard-core limit, nj > 1. 
Thus anyons can exchange positions, changing the phase of 
the total wavefunction, and show non-trivial features. 



Experimental realization 

In our experimental proposal, the key point for realizing 
anyonic statistics is to induce a hopping term with a phase 
shift which depends on the occupation of the left-hand site 
j, Fig. lb displays the basic concept. In order to distin- 
guish between different local occupational states, we require 
a non-zero on-site interaction U. For simplicity, let us con- 
sider lattice site occupations that are restricted to nj = 0, 1,2 
(higher local truncations are also possible, see Methods). 
The occupation-dependent tunneling and thus the conditional- 
hopping model equation (4) can be implemented in OLs with 
present experimental techniques. We propose to employ an 
assisted tunneling scheme, based on ideas by Jaksch and 
Zoller 19 and Juzeliunas et al. 20 . The OL is tilted, with an en- 
ergy offset A between neighbouring sites, this additional field 
gradient breaks reflection parity. Two different occupational 
states (note that the occupational state nj = is not relevant) 
in either of the two sites form in total a 4-dimensional atomic 
ground state manifold, which we propose to couple to an ex- 
cited state \e) via four external driving fields (labeled 1, 2, 
3 and 4 in Fig. lb). According to this notation, singly- and 
doubly-occupied states are coupled by fields 2 and 1 in the left 
site and by 3 and 4 in the right site, respectively. The excited 
state can be experimentally realized in at least two alternative 
ways. 

First, two spin-dependent lattices 21-23 can be employed. In 
the case of Rb 87 , one lattice for example traps atoms in the 
F = \,mp = — 1 hyperfine state, assigned to the ground state 
manifold. The excited state \e) can then be engineered as 
a vibrational state of a second lattice, trapping atoms in the 
F = l,m^ = hyperfine state. Atoms in the excited state 
would then be localized between the left and right wells of 
the F = 1 , mjr = — 1 lattice, but not necessarily in their center. 
This implementation would offer the advantage of external 
driving fields in the radio-frequency regime. Such frequencies 
could then still resolve 22 the typical energies U and A (both of 
the order of a few kHz), which is necessary to selectively cou- 
ple to the four different states in the ground state manifold. 
Second, one can use two optical lattices, and trap ground state 
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Figure 2: Density distribution in quasi-momentum space (n^) as 
a function of particle statistics 0. The shift of the density peaks 
with increasing 6 displays a characteristic quadratic behaviour. The 
fit to the trace of density maxima is depicted in white and yields 
fitting parameters k = 0.98287T, a = - 1 .051 1 /n, j3 = 0.99827T. The 
inset displays the decrease of the peak occupations with 6 (yellow 
circles) and indicates the statistically induced phase transition from 
a superfluid to a Mott-like state. Parameters: L = 30,Af = 31,(7 // = 
0.2. 

manifold atoms in the red-detuned lattice, while the excited 
state would live in the blue-detuned one. The driving fields 
necessary in this case would be then, however, typically in the 
THz frequency range, making a precise resolution of U and A 
more challenging for the experimentalist: in principle a laser 
with a linewidth cWwidth <C U, A would be needed. 
The effective tunneling rates J a b (a £ {1,2},/? G {3,4}) are 
calculated for four A- systems, laser frequencies tfy, and Rabi 
frequencies (i G {1,2,3,4}) via adiabatic elimination, see 
Methods. We emphasize that the tilt energy A disappears in 
the effective Hamiltonian after rotating out time-dependent 
phase factors: indeed this energy offset is absorbed (or re- 
leased) by the external radiation field, yielding a total Hamil- 
tonian without a tilt term (see also Jaksch and Zoller 19 on this 
issue). 

The following conditions on the effective tunneling rates J a b 
have to be satisfied in order to realize our model in equation 
(4): 

J23 = J24 = J • (5) 
J13 = Ju=Je W , (6) 

where 6 is the anyonic exchange phase. For a more detailed 
consideration of realistic energy scales and appropriate pa- 
rameter regimes, see Methods. 

Density in quasi-momentum space 

We have computed the ground state wave function for the 
conditional-hopping Bose-Hubbard model, equation (4), us- 
ing the Density Matrix Renormalization Group (DMRG) 24 ' 25 . 
In Fig. 2 we plot the quasi-momentum distribution 

(n k ) = \^ KXi ~ Xi) ^]bj) (7) 

L ij 



as a function of the statistical phase angle 6. The case 
6=0 corresponds to ordinary bosons, which for the pa- 
rameters chosen quasi-condense. The density distribution in 
quasi-momentum space is thus peaked at k = 0. Increas- 
ing 6 to non-zero values, we find that the position of the 
peaks 6 max (k) is shifted as a non-linear function of k. In- 
deed, for fillings N/L > 1 one finds a quadratic behaviour 
G max (k) = a(k - ko) 2 + /3. For fillings close to N/L = 1, 
j3 7T, ko 7T and a -> — 1 /%. For higher fillings N/L ->> 2, 
these fitting parameters are altered, however, the characteris- 
tic quadratic dependence is conserved. 

In this analysis, we find two important characteristics of con- 
ditional hopping bosons and thus anyons in ID. The quadratic 
dependence of 6 max on k contrasts ordinary magnetic fields 
(with a constant phase factor e l ° in the kinetic Hamiltonian). 
In this case, the shift of the peaks depends linearly on the 
phase angle 6. In the anyonic case, however, the growth of 
correlations with increasing may induce the characteristic 
quadratic trace, which could be directly observed in OL ex- 
periments using standard time-of-flight imaging. 

An even more important observation is that the contrast of 
the peaks (and the phase coherence) decays with increasing 
6. The peak values (n^O = 9 max )) are plotted in the inset 
of Fig. 2, in the pseudofermionic limit (6 — » n) the peak is 
almost washed out. This suggests that an increase of 6 trans- 
forms the initial quasi-condensate into a quantum state where 
all phase coherence is lost. It will become evident in the sub- 
sequent paragraphs that this quantum state will turn out to be 
a Mott-like state, with Mott plateaus emerging at fractional 
densities. We emphasize that this quantum phase transition 
is only driven by the statistical angle 0, all other parameters 
such as J/U are fixed. The loss of coherence can be under- 
stood as follows. With increasing 6 the occupation-dependent 
phase factor in equation (4) becomes more and more impor- 
tant: Hopping processes connecting sites with different oc- 
cupations will contribute different phases and will cancel out 
in the kinetic Hamiltonian due to incoherent superpositions. 
This destructive interference effect is amplified by an increas- 
ing 6 and induces the localization of particles, yielding an in- 
sulating phase. We emphasize that the present analysis of the 
density distribution in momentum space refers to the bosonic 
particles only, (tik) = {b\bk). Namely, Fig. 2 represents what 
would be observed in the experiment. However, while the 
mapping in equation (2) establishes a 1-1 correspondence be- 
tween the number operators in real space, (a^aj) = (b^bj), 
the density distributions in momentum space may differ sig- 
nificantly. In this sense, a study of % = (tfj^fc) an d the su- 
perfluid order parameter in the original anyonic model (3) 
would be very interesting, but it is beyond the scope of this 
paper. A study of density distributions in momentum and real 
space was recently presented for the particular case of hard- 
core anyons 26 . 



Phase diagram 

Next, we present the phase diagram for conditional- 
hopping bosons in the (fl^J/U) -plane. Just as in the case of 
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ordinary bosons (0 = 0) - the celebrated phase diagram of the 
Bose-Hubbard model 9 - the anyonic version will display insu- 
lating and superfluid regions. However, we find that the size 
of the insulating regions (Mott lobes) grows with increasing 
statistical angle 6 - a fact that will be central in our proposal 
to design phase transitions by tuning the particle statistics. 
The phase diagram is calculated as follows. We start with a 
unit filling ground state \N = L) and compare the energies with 
the ground states of \N = L± 1). This yields the gap energies 
AE ± (L), corresponding to the energy required to add or sub- 
tract one particle, respectively. These gap energies were cal- 
culated using DMRG for system sizes L = 15,30,60,90 and 
120. From the finite-size scaling we extrapolate the infinite- 
size values AE ± , which are plotted in Fig. 3 in the ( jix, //£/)- 
plane. Note that a Mott insulator is associated with a non-zero 
gap £ = AE+ — AE~, while for the superfluid phase £ = 
in the thermodynamic limit. For the bosonic case = we 
recover the first Mott lobe as in the work by Kiihner and 
Monien 27 . Turning on the statistical angle and transmuting 
the bosons into anyons, we observe an expansion of the Mott- 
like insulating phase in both dimensions jl and J/U. The 
phase transition points (U /J) C Ht (defined by the cusps of the 
Mott lobes) are plotted in the inset of Fig. 3 as a function 
of 6. In the pseudofermion limit 6 —> 71, the Mott- like phase 
seems to extend to very large values of J/U. In contrast, or- 
dinary bosons would form a superfluid state in this parameter 
regime. The expansion of the Mott lobes with 6 is also ob- 
served in mean-field calculations for our model, equation (4). 
The mean field solution produces very interesting patterns for 
the transition lines, as shown in Fig. 5. For details of this cal- 
culation, see Methods. 

We note that the phase diagrams for conditional-hopping 
bosons and for anyons are the same. Due to the unitarity of 
the mapping (2), the two models (3) and (4) are isospectral, 
and thus feature the same energy gaps and phase diagrams. 



Discussion 

We envision the following procedure to demonstrate the 
first statistically induced phase transition: We fix the param- 
eters at J/U = 0.5, N/L = 1 and \i ~ 0. We start with a phase 
detuning = between the external driving fields (see Meth- 
ods), and thus realize a superfluid bosonic gas as the ground 
state. The phase angle 6 is now continuously increased, lead- 
ing to anyonization of the gas and growth of the Mott-like 
phase. At a critical value 6 C the phase border will surpass the 
fixed point in parameter space, which will be then located in- 
side the Mott phase. The critical angle can be estimated from 
the phase diagram to be in the range 6 C G [;r/2,3;r/4]. For 
laser detunings beyond this critical angle, 6 > C , the gas will 
be in an insulating, Mott-like state. In this way, the varia- 
tion of the particle statistics in our proposal directly realizes a 
novel superfluid-to-Mott quantum phase transition. 

An intriguing aspect of the Mott-like state emerges when 
a harmonic trapping potential is added to the system, H% r = 
H b + V l)/2 — i) 2 rti. This simulates the experimental 
conditions also to a more realistic extent. In a harmonic trap, 
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Figure 3 : Phase diagram of the anyonized gas. The Mott lobe (cor- 
responding to (n) = 1) expands with increasing statistical angle 6 in 
both directions in the (fJL,J/U) -plane. This demonstrates the novel 
possibility to induce a quantum phase transition from the superfluid 
into the insulating, Mott-like phase, by simply changing the particle 
statistics. The energies AE ± are symbolized by circles (crosses), re- 
spectively. The critical phase transition points (U /J) C rit are shown 
on the inset, as a function of 6. In the extreme limit — > 71, (U /J) cr it 
tends to zero, i.e. pseudofermions seem to be always in the insulating 
phase. 



ordinary Mott insulators form real- space density distributions 
(rii) that resemble "wedding cakes" 9 . Due to the broken trans- 
lational invariance (induced by the trap), the chemical poten- 
tial in the local density approximation now is a function of 
the lattice sites. The distribution (n/) thus exhibits plateaus 
at integer densities, for ordinary bosons in a Mott-insulating 
state. 

We have computed the density distribution for a system 
with additional harmonic confinement, at the same fixed pa- 
rameter point as in the procedure outlined above. The result is 
plotted in Fig. 4. For low statistical phase angles 6 (and thus 
a superfluid state), the density distribution displays a smooth, 
quadratic profile centered around the minimum of the trap. 
Beyond the critical phase angle, 6 > C , Mott-like plateaus 
appear in (n/). Surprisingly, in addition to the integer density 
Mott plateau at n = 1, a new plateau emerges at the fractional 
value n ~ ^ . Fractional Mott plateaus persist also in other pa- 
rameter regimes, when J/U and the trapping potential V are 
varied, and seem to form a universal feature of ID anyonic 
gases at large exchange phases 0, which will be subject of a 
future study. 

Recent progress on the experimental side have made direct 
imaging of density distributions possible, using "quantum gas 
microscopy" 28 ' 29 . A few weeks ago, the plateau structure of 
the Mott insulator was directly observed for the first time 30 at 
the single atom level. This new technique opens up the possi- 
bility to directly demonstrate fractional Mott plateaus and to 
image statistically induced phase transitions in situ. 

In summary, we have shown how to realize fractional statis- 
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Figure 4: Fractional Mott plateaus. The density distributions (ni) 
are plotted for increasing statistical angle 6. In each panel the density 
distribution for ordinary bosons (6=0, blue line) is drawn for refer- 
ence. The distribution at high statistical angles (6 = ^,7T, marked 
black and green) exhibit a Mott-like plateau at the fractional value 
n ~ \ (marked by arrows). Parameters: N = L = 30, J/U = 0.5, 
V/U = 0.0l. 



tics in ID optical lattices, using bosons in a realistic experi- 
mental environment. The experiment we propose features the 
full control and tuneability of the particles' exchange statis- 
tics - paving the way to the first statistically induced quantum 
phase transition. 



Methods 

Realizing conditional-hopping bosons 

In this section we discuss how to realize the conditional- 
hopping Hamiltonian equation (4), using four different and 
independent A- transitions 31 . In general, we assume a deep 
optical lattice potential, giving rise to a negligible bare kinetic 
tunneling amplitude J^ n . Let us focus for a moment on one 
A-scheme, where two ground state levels \a) and \b) are cou- 
pled through a Raman process via an excited state \e). In our 
case, | a) and \b) correspond to the wavefunctions of atoms 
at distance d localized in neighbouring sites of the tilted lat- 
tice V(x) (d being the lattice constant), while \e) experiences 
another potential V f (x) (for a brief discussion of two realistic 
experimental possibilites to realize V'(x), see the main text). 
The levels have energies E\ = ftCQi, i = a,b,e and the transition 
between a(b) and e is driven by an external radiation field with 
frequency CO e — CO a ^ — 8, with detuning 8. Note that in our 
scheme the energies E[ depend on A and U . The Hamiltonian 
for the three-level A-system reads 



H= £ h(O i \i)(i\ + -(Ya\e)(a\ + Y b \e)(b\+h.c.) (8) 

i=a,b,e 



where y a{b) = ^ w ^ ) e" I ' (fflfe_fi)fl W" 5)/ . Here, the quantity 
is the Rabi frequency for the transition a(b) —> e with the 



atom centered in the same position. However, since ground 
and excited states feel different lattices, the x components of 
the Wannier functions w(x) are slightly displaced. Thus, the 
off-diagonal elements in equation (8) contain the integrals 19 

W e a = e ikaXa f w* e (x + x e )e ikaX w a (x)dx, 

w e = e ik b (x a +d) J w *( x+Xe )^ Wb ( x + d)dx, (9) 

where k a ^ is the x-component of the momentum of the driv- 
ing radiation field, with modulus \k a (p) | = (co e — ® a (b) ~$)l c - 
x a is the atom position in the left well, while x e refers to the 
position in the excited state. Note that the integrals in (9) are 
in general non-zero, as the two Wannier functions belong to 
different lattices. The quantities y a (b) are complex numbers, 
whose modulus and phase can be freely tuned by choosing the 
appropriate intensity, polarization and direction of the driv- 
ing fields. For sufficiently large detunings 5 > \y a (b)V the 
level \e) can be adiabatically eliminated and in the rotating 
wave approximation the effective Hamiltonian in the subspace 
{\a), \b)} reads 



#eff : 



ft 



( 



\fa\ 2 



Ta% 



45 V Tbfa \?b\ 2 



(10) 



where f is the non-rotating version of 7, i.e. without the time- 
dependent phase factors. In order to realize equation (4), we 
suggest to employ four driving fields (see labels in Fig. lb) 
with different frequencies in order to avoid interferences. This 
situation would correspond to a maximum of two atoms per 
site. If the local density truncation were set to a higher num- 
ber, more external driving fields would become necessary. As 
all fields can be tuned independently from each other, this 
poses no problems besides potential budgetary considerations. 
The couplings /13, Ju, J23 and J24, between the four different 
levels are then obtained in terms of the effective Rabi frequen- 
cies, J ab = %%/28. 

Our aim is to satisfy the conditions set in equations (5-6) in 
order to engineer a state-dependent phase factor. This im- 
plies 72 = f\e ld , which can be achieved by the free tuneabil- 
ity of each driving fields 's frequency, intensity, polarization 
and direction. Furthermore this choice of parameters implies 
\fa\ = \%\, i-e. the diagonal elements of the effective Hamil- 
tonian are now equivalent. Thus, the tilt energy A has van- 
ished via adiabatic elimination, and consequently also does 
not appear in equation (4). Indeed this offset energy between 
neighbouring sites is compensated by the external radiation 
field. As the assisted tunneling proposed in this Article is the 
only mechanism for hopping in the lattice, unwanted effects 
like Bloch oscillations do not appear in our system (the bare 
kinetic tunneling amplitude J^ n is assumed to be negligible 
compared to all energy scales discussed here). 
In summary, the parameters discussed here have to satisfy the 
following conditions: 

First, Siinewidth <A,t/, so that the external driving fields re- 
solve the different levels of the ground state manifold. Sec- 
ond, large detunings 8 > \Y a (b) \ are required for a short-lived 
excited state and the validity of the adiabatic elimination. 
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Third, A and U can be in the same frequency regime (a few 
kHz), but should not be identical (their difference should 
be ^> <5iinewidth)- As an example, A ~ 2 kHz, U ~ 3 kHz, 
\J ab \ = J ~ 5 kHz and |y fl (^)| — 20 kHz would be sufficient 
if the linewidth of the radiation field were 5ii new idth — 50 Hz, 
which is a realistic assumption for typical radio-frequency 
driving fields (see e.g. the works by Campbell et al 22 and 
McKay and DeMarco 23 ). 



Fractional Jordan-Wigner mapping 

In the following we prove that anyons are isomorphic to 
bosons in ID. 

In particular we prove that the operators a, as defined in the 
non-local mapping, equation (2), indeed obey the anyonic 
commutation relations of equation (1), provided that the par- 
ticles of type b are bosons. 

For the case i < j we wish to rewrite products of anyonic 
operators in terms of the bosonic ones: 

a/at = bie- WLi < k <J nk b] 
1 j 1 j 

= e- ie ^< k <J nk bib]e- ieni , 
f{p)a)oi = e- ie ^< k <j nk e- ieni f(e)b)bi 

= e -i6Isi<k<jnk e -iQ(ni+\)y\fo^ ^\) 

Here we have defined f(0) = e Ws ^-^ and used that f(0) = 
e -i0 smce i < j was assumed. We can now evaluate the LHS 
of equation (1): 

aia)-f(e)a)ai = e - ie ^< k <J nk (bib)e- Wni - e - ie ^%)bi) 

= e - ie ^< k <j nk e- w ^ l) [bi,b]} 

= (12) 

Thus the anyonic commutation relations have been proven for 
the case i < j. The proof for the case i > j is very similar. 
For the case i = j one just has to note that a] at = b\b[ and 
f(6) = 1. Note that the resulting conditional-hopping bosonic 
Hamiltonian, equation (4), resembles the exactly solvable 
model of q-bosons 32 . However, this model is not equivalent 
to our model. 



Mean-field calculation 

The conditional-hopping bosonic Hamiltonian is given by 
equation (4). In addition, we now include also a chemical 
potential term. Expressing energies in units of U (we fix U = 
1), we have 



6=t:/4 



6=n/2 



H = L \ u Mj ~ l ) - ^ n J - J ( c ] b j+i+ b )+i c j) 



(13) 
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Figure 5: Mean-field solution for the Mott-superfluid transition. 

The transition lines are displayed for several values of 6. For com- 
parison we display also the data for the bosonic case 6 = (blue 
curve). 



In absence of hopping / = 0, all the sites are independent 
and the ground state is of the Gutzwiller type 



l*Po> = |v>* 



v= Vv! 



where v = N/L is the filling factor. The local energies are 
3 1 



e(v) = iv(v — 1) — /IV, while the local gaps for adding or 



subtracting one boson are, respectively 

e(v + l)-e(v) = v-ju, e(v-l)-e(v) = - (v-l)+/i, 



where for convenience we have defined cj = e lQn ibj. 



So, the ground state in every site has v particles in the interval 

^li v) < \i < ^l| v) , with jui v) = v - 1 and ^l| v) = V. The gap 
in the whole system at a given number of particles is given by 
A = 1 , obtained by removing an atom at some site an putting 
it in another one already occupied. 

Here the mean field (MF) is obtained by decoupling the 
hopping term as Cjbj+i « — a|ai + (X^bj+i + (XiCj, where the 
order parameters are OC\ = (bj) and Ct2 = (cj). Accordingly, 
the Hamiltonian (13) in MF becomes 



with 

Hj = l -nj{rij - 1) - \inj-jYj^ib) + o%bj + a x c] + a*Cj). 
1 j 

The two parameters (X\ and a 2 are not completely indepen- 
dent as as they are both vanishing or non vanishing at the 
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same time. Of course there is a trivial solution correspond- 
ing to OL\ = 0, / = 1 , 2, which corresponds to the Mott insu- 
lating phase. The occurrence of a/ ^ signals the instability 
towards superfluid correlations. On inhomogeneous lattices 
a further situation can in principle occur, where OL\ ^ only 
on a fraction of the lattice sites. The self-consistent relation 
defines a map a/ = A/// ay, obtained by linearizing about the 
solution a/ = 0. The instability of the trivial solution sets in 
when the maximal eigenvalue of A is greater than 1 . Close to 
the trivial point, it holds \a\\ <C 1, hence the kinetic term can 
be treated perturbatively. Up to first perturbative order, the 
wavefunction can be written as \\jf) = -f Ji// 1 )), where 

|yA(°)) = |v)and 



V e(v)-e(v') |V) 

r V^(a| + a i^ 0(v_1) 



/i-v + 1 

Vv + l(a 2 + ai^ 0v ) 
v-/x 



v-1) 
v + 1) 



Hence using the self-consistency relations GC\ = {\\r\bj\\\r) 
and a 2 = (\jf\c , 



ai/J 
a 2 /J 



v(a 2 + aie w ( y -^) (v + l)(a 2 + ct\e m 



/i-V + 1 V — \i 

v(a 2 ^' ( v - 1 ) + ai) (v + l)(a 2 ^' 0v + ai) 



The matrix A is then 



with 



A + (^' e -l)5 



( M -[ M ])([ M ]- M + 1)' 



5: 



since every lobe is labelled by v = \\\\ + 1. The eigenvalues 
of A are given by 



/(e) +/(-e) ± + (/( ) _/(_e))2 



The condition max{A + ,A_} = 1 signals the onset of insta- 
bility of the trivial solution and hence establishes the critical 
coupling J cr it along the Mott- superfluid transition line. In Fig. 
5 the phase diagrams for several values of 6 are shown. The 
Mott lobes expand for non-zero statistical angles 0, a fact cen- 
tral for designing statistically induced phase transitions. 
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